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Abstract 



We consider the effects connected with the detailed radiative transfer during the epoch of cosmological recombination on the ionization 
history of our Universe. We focus on the escape of photons from the hydrogen Lyman a resonance at redshifts 600 < z < 2000, one 
of two key mechanisms defining the rate of cosmological recombination. We approach this problem within the standard formulation, 
and corrections due to two-photon interactions are deferred to another paper. 

As a main result we show here that within a non-stationary approach to the escape problem, the resulting correction in the free electron 
fraction, N^, is about ~ 1.6 - 1.8% in the redshift range 800 < z < 1200. Therefore the discussed process results in one of the largest 
modifications to the ionization history close to the maximum of Thomson-visibility function at z ~ 1 100 considered so far. 
We prove our results both numerically and analytically, deriving the escape probability, and considering both Lyman a line emission 
and line absorption in a way different from the Sobolev approximation. In particular, we give a detailed derivation of the Sobolev 
escape probability during hydrogen recombination, and explain the underlying assumptions. We then discuss the escape of photons 
for the case of coherent scattering in the lab frame, solving this problem analytically in the quasi-stationary approximation and also 
in the time-dependent case. We show here that during hydrogen recombination the Sobolev approximation for the escape probability 
is not valid at the level of AP/P ~ 5 - 10%. This is because during recombination the ionization degree changes significantly over 
a characteristic time Az/z ~ 10%, so that at percent level accuracy the photon distribution is not evolving along a sequence of 
quasi-stationary stages. Non-stationary corrections increase the effective escape by AP/P — [-6.4% at z ~ 1490, and decrease it by 

AP/P 7.6% close to the maximum of the Thomson-visibility function. We also demonstrate the crucial role of line emission and 

absorption in distant wings (hundreds and thousands of Doppler widths from the resonance) for this effect, and argue that the final 
answer probably can only be given within a more rigorous formulation of the problem using a two- or multi-photon description. 

Key words. Cosmic Microwave Background: cosmological recombination, temperature anisotropies, radiative transfer, theory 



1. Introduction 

The extraordinary advances in observations of the Cosmic 
Microwave Background (C MB) tempe r ature a nd polarization 
angul ar anisotropies (e.g. iPage et aT] 120061: iHinshaw et aLl 
|2006|) and the prospects with the Planck SurveyoiQ have 
motivated several groups to re-examine the problem of 
cosm ol ogical recombina tion (e.g. see ISunvaev & Chlubal 
120081: iFendt et all l2008l for detailed overview), includ- 
ing subtle physical pr ocesse s during hydrogen (e.g. see 
Dubrovich & Grachev 120051: |Chluba & Sunyaev 2006b; 
Kholucenko & Ivanchik 2006; Rubino -Martfn et al.1 120061: 
iChluba& Sunyaey 2007; Hirat a 2008) an d h elium recombina- 
tion (e. g. see ISwitzer & Hiratal |2008all5: iHirata & Switzeij 
_lKholuDenko et al.' '2007"; 'Wong & Scotd l2007t 
Rubino-M artfn et al. 2008; Kholupenko et al. 200^. It has 
that percent level corrections to the ionization 
which should be taken into account for future 
parameters using CMB data 



been argued 
history exist, 

determinations of cosmological 
obtained with the Planck Surveyor 

In this paper we investigate the validity of one of the key 
simplifications used for computations of the hydrogen recombi- 
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nation history within existing multi-level recombination codes: 
the Sobolev approximation for the escape of Lyman a photons 
from the center of the resonance. With this approximation it is 
possible to separate the problem of the evolution of the photon 
field and the populations of the hydrogen atom. 

Originally the Sobolev approximation was developed in 
order to describe the escape of photons fro m finite expand - 
ing envelopes of planetary nebulae and stars (ISobolevlll960l) . 
but it has been shown that even for cosmological applica- 
tions i.e. infinite slowly expand i ng media, it is very use- 
ful (iQrache v & DubroyichI 1 1 9911 ; iHummer & Rvbickil 119921; 
Rybicki & deirAntonio!ll994|)~ It gives the same answer as less 
sophisticated methods, based on simple solutions of the integral 
or differential equations of radiative transfer, which were used to 
solve the cosmolo gical hydrogen reco mbination pro blem in the 
1960s dVarshal ovich & SvunvaevI [T968: Zeldovich et al.l [19681; 
Peebles 1968). Both for the Sobolev approximation and these 
simpler derivations the main assumptions are: (i) the properties 
of the medium (e.g. ionization degree, density, expansion rate) 
do not change much while the photons interact strongly with the 
Lyman a resonance and (ii) each scattering leads to a complete 
redistribution of photons over the whole line profile. 

Due to assumption (i) it is possible to approximate 
the evolution of the photon distribution as quasi-stationary, 
which for conditions in our U niverse seems to be reasonable 
(iRvbicki & delFAntoniol 19941) . However, up to what level of ac- 
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curacy remains a difficult question and deserves further investi- 
gations. On the other hand, assumption (ii) is much less justified, 
since complete redistribution requires some process that destroys 
the coherence in the resonance scattering event. This is usually 
done by collisional processes, which for t he conditions in ou r 
Universe are extremely inefficient (e.g. see lChluba et alj|2007l) . 
We will demonstrate here that for present day experimental re- 
quirements, i.e. sub-percent level accuracy in the theoretical pre- 
dictions of the CM B power spectra at large multipoles I (e.g. see 
ISehak et"ani2003h . both approximations become insufficient. 

In order to understand this problem, it is important that the 
ionization degree during cosmological hydrogen recombination 
changes with characteristic time Az/z ~ 10%. Also, it is clear 
that photons, which are released in the distant wings of the 
Lyman a linfl can in principle travel, scatter, and redshift for 
a very long time until being reabsorbed. Here it is important to 
distinguish between line scattering events, and line emission and 
absorption processes. The former only lead to a redistribution 
of photons over frequency, but no net change in the ionization 
degree, while the latter (which for example are connected with 
direct transitions of electrons between the continuum and the 
2p state) are able to change the number of Lyman a photons, 
and hence the ionization degree. Note that during hydrogen re- 
combination, line absorption occurs with much lower probabil- 
ity (~ 10""^ - 10"^) than line scattering, so that photons only die 
or disappear effectively rather close to the line center (within a 
few ten to hundred Doppler widths from the resonance), while in 
the distant wings they mainly scatter. In addition, every photon 
that was absorbed (or died) will be replaced by a new photon in a 
line emission event after a very short time. The profile of this line 
emission is usually described by a Voigt profile, so that the com- 
bination of hne absorption followed by a line emission appears 
to lead to a complete redistribution of photons over the whole 
Lyman a line profile. However, from the microscopic point of 
view this is not a scattering event0. 

As explained in iRvbicki & delF Antonio! (1 19941) . in the ex- 
panding Universe the redistribution of photons due to Lyman 
a resonance scatterin g is more accurately described by so-call 
type-II redistribution (lHummerl ll962) rather than by complete 
redistribution. In the former case the photon scatters coherently 
in the rest-frame of the atom, so that the changes in the ener- 
gies of the photon after the scattering event are related to the 
motion of the atom. Studying this type of redistribution pro- 
cess in de tail is beyond the scope o f this paper, but our com- 
putations (IChluba & SunvaevI l2009h. in very good agre ement 
with earlier works (e.g. see^ Rvbicki & deirAntonioll994l) . show 
that in a time-dependent formulation of the problem, including 
Doppler broadening, atomic recoil and stimulated emissio 
actual solution for the scattered photon distribution is very close 
to the one in the case of no redistribution, or equivalently no 
line scattering. Here we show in addition that the assumption 
of complete redistribution leads to several unphysical conclu- 
sions, both in the quasi-stationary approximation and a time- 
dependent approach. This is due to the very peculiar conditions 
in our Universe, where collisional processes are not important. 



^ At redshift z = 1100 a thousand Doppler width corresponds to 
Ay/v ~ 2%, a distance from the line center that can be passed by red- 
shifting in Az/z ~ 2%. 

^ ISwitzer & Hiratal l l2008ah also make this distinction using the ter- 
mini of incoherent processes and coherent scattering. 

Within a Fokker-Plan ck approach the atomic recoil effect was first 
included bv lBaskd il978l) . while the eff ect of stimulate d emission was 
only taken into account very recently bv lRvbickil ( l2006h . 



and in particular where due to the low Hubble expansion rate, 
the Sobolev optical depth reaches extreme values of ~ 10^-10^ 
during recombination. 

We therefore investigate the evolution of the photon field in 
the no-scattering approximation and show that time-dependent 
corrections to the effective escape probability are important at 
the level of ~ 5% - 10% (see Sect. [Ml Fig. Hand Fig.©. As 
mentioned above, this is due to the fact that in the distant wings 
of the Lyman a resonance photons mainly scatter, but do not 
disappear. This renders it important to include changes in the 
ionization degree and photon emission rate during the evolution 
of the photon field in the computations, implying that the quasi- 
stationary approximation becomes inaccurate. Both changes in 
the absorption optical depth and the effective emission rate can- 
not be neglected. The coiTesponding time-dependent changes in 
the free electron fraction, which are important for the Thomson 
visibility function and in computations of the CMB power spec- 
tra, reach the level of ~ 1.6 - 1.8% in the redshift range 800 < 
z < 1200 (see Sect. |4]and Fig. [T2l i. and therefore are about 2 
tim es as large as tho se due to atomic recoil, recently studied by 
Grac hev & Dubrovi ch (2008). Taking the time-dependent cor- 
rection investigated here into account will therefore be very im- 
portant for the analysis of future CMB data from the Planck 
Surveyor. 

We also briefly discuss another aspect of the Sobolev ap- 
proximation, which is connected to the shape of the Lyman a 
line profile (see Sect. 13.5b . In the Sobolev approximation there 
is no direct dependence of the Sobolev escape probability on the 
shape of the line emission, absorption, and scattering profiles, as 
long as all are identical. Our derivation also clearly shows this 
point (cf. Sect. 13.21 and 13. 5l l. Therefore, in principle it does not 
matter if the profile is assumed to be a Lorentzian, a Voigt pro- 
file, a pure Doppler profile, or a (5 function. It also turns out that 
in the no line scattering approximation this is true, as long as the 
line emission and absorption profiles are identical, and the evo- 
lution of the photon distribution is quasi-stationary (cf. Sect. 13.31 
and[331l. 

However, if one includes the deviations from quasi- 
stationarity, then the result does depend in detail on the Lyman a 
profile, even if the line emission and absorption profiles still are 
the same. For example, in the case of a pure Doppler profile (very 
naiTow), the problem of the Lyman a photon escape from the 
resonance due to the expansion of the Universe would be prac- 
tically quasi-stationary, and the Sobolev approximation should 
be applicable. This is because the number of photons emitted 
and absorbed in the optically thin region of the Lyman a line is 
exponentially small, and all the transfer is happening inside the 
Doppler core, coiTesponding to Av/v ~ few x 10"^. 

On the other hand in the real problem, Lyman a emission and 
absorption also occurs in the distant Lorentz wings (at hundreds 
and thousands of Doppler widths) of the resonance. As we show 
here, at a percent level the number of these photons is very im- 
portant for the value of the effective escape probability (e.g. see 
Fig. [Tol l. This shows that it is crucial to understand the profiles 
(or cross-sections) of the considered processes in more detail, 
and for this probably a formulation in the two- or multi-photon 
picture will become necessary. Also in principle it should be pos- 
sible to distinguish between different redistribution processes for 
the line scattering event, by measuring the shape and position of 
the residual, present day CMB Lyman a distortion. 

It is extremely impressive that the standard estimates of the 
Lyman a escape probability, which were used in the first papers 
on cosmological recombination, and the Sobolev approximation 
give such precise (better than 5-10%) answers, even though they 
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are based on two incorrect assumptions as mentioned above. It is 
well known that the principal difference (from a physical point 
of view) between the cases of partial and complete redistribu- 
tion does not influence the fin al result very much in the majority 
of astrophysical applications (llvanovlll973l) . However, the enor- 
mous requirements of accuracy of theoretical estimates in the era 
of precise cosmology change the situation, and force us to search 
for percent level corrections to the escape of Ly a photons from 
resonance during the epoch of cosmological recombination. 

2. Transfer equations for the photon field 

In this Section we provide the transfer equation describing the 
evolution of the photon distribution in the vicinity of the Lyman 
a resonance. We include the effect of line emission and line ab- 
sorption in the expanding Universe for the cases of coherent line 
scattering in the lab frame, and complete redistribution. Here 
we envision all processes as 1 -H 1 ph oton processes, as in the 
Seaton-cascade description ( ISeatonl ll959). but leave the treat- 
ment of correction due to two-photon interactions for a future 
paper Also the effects of partial frequency redistribution will be 
discussed in separate paper. In Sect. |2.3| and |2.4| we give the time- 
dependent solutions of these equations. We will use these results 
in Sect. |3]to deduce the Lyman a escape probability, which then 
can be utilized to estimate the corrections to the cosmological 
ionization history. 

2.1. General kinetic equation for the photon field 

To follow the evolution of the photon field in the expanding 
Universe we start with the kinetic equation for the function Ny - 
ly/hv, where ly is the physical s pecific intensity of the isotropi c, 
ambient radiation field (e.g. see lRybicki & dell' Antoniol 19941) : 



dNy 



dt 



+ 2HNy-Hv 



dNy 
dv 



■ C[Ny] . 



(1) 



Here H(z) is the Hubble parameter as a function of redshift z 
and C[Ny] is the collision term, which describes the emission, 
absorption and frequency redistribution processes. 

In order to simplify the left hand side of the equation ^ we 
transform to the frequency variable x - v/(l + z), so that 



dx 
dv 



(1+z) 



Inserting this into Eq. (HJ yields 



ON, 



dt 



+ 3HN^-Hx — i- 
ox 



^{l+Z)C[Ny]. 



To obtain dN^/ dt\y one can use the total differential of Nx 

ON, 



dN, 



dt + 



dx 



dx 



which with dxj dt\y - xH then gives 



dNx 
dt 



dt 



+ Hx 



dx 



Inserting this into Eq. ([3]) one finds 

dN 



dt 



+ 3HNx 



= il+z)C[Ny]. 



(2) 



(3) 



(4) 



(5) 



(6) 



Here the redshifting term was absorbed due to the choice of the 
frequency variable. 

The term 3HNx can be eliminated using the substitution 
ff^ - NJ(l + z)^ = Ny/(l + zf, so that Eq. O takes the form 



1 dN, 
c ~dt' 



C[Ny] 
(l+z)2 



(7) 



One can easily verify that in the absence of physical interactions 
(C[Ny] = 0), in spite of the Hub ble expansion, a Pla nckian spec- 
trum is not modified (e.g. see iPadmanabhanll2002l) . so that it is 
always possible to directly write 



1 dN, 

c dt 



1 dANx 
c 



dt 



(8) 



where A^V^ - Nx - N^^ is the corresponding deviation of the 
spectrum from a blackbody, which in our coordinates reads 



Nl' = ^ 



2 ghx/kTo _ I ■ 



(9) 



Here Tq = 2.725 K is the CMB temperature today 
(iFixsen & Mathedl2002h . 



2.2. line emission and line absorption 

Although for conditions in the Universe during cosmological re- 
combinatior0 the resonant scattering optical depth close to the 
Lyman a line center exceeds unity by several orders of magni- 
tude, only real line emission and absorption lead to a net change 
of the photon number If we consider an electron in the ground 
state of hydrogen which after the absorption of a photon (say 
close to the Lyman a resonance) is excited to the 2p state, then 
there are two routes out of this level: (i) it can directly decay back 
to the ground state, re-emitting a photon with (slightly) changed 
frequency, depending on the considered redistribution process, 
or (ii) it can be further excited to the continuum or higher shells 
(n > 2) by the subsequent absorption of a blackbody photon 
from the CMB. Only in case (ii) does the number of Lyman a 
photons really change, while for (i) the photon simply was scat- 
tered. 

To describe this aspect of the problem, we use the death 
probability or single scattering albedo, p^, which specifies what 
fraction of photons that interact with a hydrogen atom in the 
Is state, will really disappear from the photon distribution. The 
scattering probability, Psc = 1 - Pd, will then give the fraction 
of photons that reappear at a different frequency, and hence only 
underwent a scattering rather than a real line absorption. 



2.2.1 . Death probability or single scattering albedo 

Including all possible ways in and out of the 2p level, the net 
change in the number density of electrons in the 2p level can be 
written as 



djV2p 
dt 



+ 3HN2r, = 



dt 



hy-a 
2p 



+ ^2p-^2p^2p' 



(10) 



' Electron and proton collisions are negligible in comparison with 
radiative processes, like photorecombin ation or photoionization, and 
bound-bound dipole transitions (e.g. see iChluba & Sunvaevll2008bh . 
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dA^ iLy-Q- 

where —^\ denotes the contribution from the Lyman a tran- 

l2p 



dl 



sition, which we will specify below (see Sect. 13.1b . and 



^2p ='R2pc + V — A/2pnP'(vap) 

.?2p 



/>2p 



(11a) 



/?+p = A^e A^p Rc2p + ^'2p [ 1 + nP'(vap)] . (lib) 

/>2p 

Here A^e and A^p are the electron and proton number densities, 
and nP' = 1 /(e'"'/*^^'' - 1) is the blackbody photon occupation 
number, where Ty = To ( 1 +z) is the CMB temperature at redshift 
z. Furthermore, A,2p denotes the spontaneous dipole transition 
rate from level ; to the 2p state, v,2p the corresponding transition 
frequency, and gi the statistical weight of level /. R^p and /?2pc 
are the photorecombination and photoionization coefficients of 
the 2p state, which are computed assuming that the ambient ra- 
diation field is Planckian. Since except for the Lyman series and 
the 2s- Is two-photon transition, the emission of photons dur- 
ing cosmological hydrogen recombination only produces tiny 
deviations of the photon dis tribution from a blackbody spec- 
trum By (e.g. see IChluba & SunvaevI l20()6ah . this approxima- 
tion is possible. Similarly, we have neglected the spectral dis- 
tortion in the terms due to transitions from and to higher levels 
(A,-2p riy « A/2p Hy'), so that R^^ becomes completely independent 
of the solution for the photon field. However, note that /?2p still 
depends on the solution for the populations, A^,, of the excited 
levels, and the electron and proton number density. 

Omitting electron and proton collisions, the total probability 
for Lyman a emission /?em is therefore given by 



(12) 



A21 +R 



2p 



where A21 = 6.27 x 10** s"' is the spontaneous 2p-ls transition 
rate. The corresponding probability for the death of photons, i.e. 
removal of Lyman a photons or return of 2p electrons to the 
continuum or higher levels, is then given by = 1 - Pem- 

Note that in Eq. ( fT2l l we directly neglected the effects of 
stimulated emission. This approximation is well justified, since 
close to the Lyman a transition the photon occupation number 

«v = TT-T 1 at all relevant redshifts. 

' 2hv-' J—. 

In Fig. \l\ we show the death probability, p^, as a function 
of redshift considering a 2, 3 and 10 shell hydrogen atom. It is 
clear that the largest contribution to the death probability comes 
from the third shell, and cases with « > 3 are practically indistin- 
guishable. This is because during cosmological hydrogen recom- 
bination /?2pc « i[A3s2p+5A3d2p]n'''(v3s2p), andsincenP'(v3s2p), 
for « > 3, is exponentially larger than n'''(v„s2p), so that also 

[A3s2p +5 A3d2p] n'''(v3s2p) » [A„s2p +5 A„d2p] nP'(v„s2p)- This fact 
implies that for a consistent investigation of the Lyman a escape 
problem, one should include at least 3 shells in the computations. 

2.2.2. Line emission profile 

The form of the emission profile for the Lyman a line (under the 
assumption of complete redistribution) is known from quantum- 
mechanical considerations. Including the thermal motion of the 
hydrogen atoms it is usually described using the so-called Voigt 
profile: 



(p(v) 



ji^l^ Avi 



-f 



e-' At _ ^(v) 



(13a) 




2000 



Figure 1. Different probabilities for the cosmological hydrogen 
recombination problem as a function of redshift. The death prob- 
abilities, pd, for a 2, 3 (dashed lines), and 10 shell hydrogen atom 
(dotted line) are shown. The death probability for the 3 shell case 
already practically coincides with the death probability for the 
10 shell case. The solid lines show different escape probabilities. 
Ps = [1 - e^^^]/Ts denotes the normal Sobolev escape probabil- 
ity, while -Pd = [1 ~ e '^''1/Td, where = Pd^s- To compute ts 
we used the Recfast solution for A^is. 

where for the H i Lyman a transition the Voigt parameter, a, the 
Doppler width of the line due to the thermal motion of the hy- 
drogen atoms, Avd, and the variable xu are defined by 



A21 





47rAvD 


Avd _ 




2kBT, 


V21 


i 






V-V21 

Avd 



8.61 X 10" 



2.35 X 10" 



(1+z) 



-1/2 



1100 a- 
(l+z) 



1100;t 



1/2 



(13b) 

(13c) 
(13d) 



Here V21 ~ 2.47 x 10'^ Hz is the Lyman a transition frequency 
and;t' = Ty/T,,. 

The Voigt profile is normalized such that J^^ ^ dvdQ = 1 
and it has the well known limiting cases 



<P(y) 



for |xd| > 10 
for Xd ~ 



(14) 



with TV = e°^Erfc(fl) ^ 1 - 2a/ a/tt + a^ for a 1 and where 
Erfc(x) denotes the complementary error function. 

In addition, on the red side of the resonance one can approx- 
imate the integral = (p(y') dv' by 

1-1/2 



Swings 



2.73 X 10 



-6 



(1+z) 



1100;r 



Xd 



-100 



(15) 



as long as -V2i/Avd xd < -10. This formula shows that 
only a very small fraction of photons is directly emitted in the 
distant wings. However it is also known that escape of photons 
from the Doppler core is extremely strongly suppressed. As a 
result the emission of the photons in the distant wings should be 
considered carefully. 
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2.2.3. Line emission term 

With the definitions given above, the term for real emission of 
photons due the addition of fresh electrons to the 2p state can be 
written as 



1 dNy 
c dt 



Ly-a 



(16) 



The emission probability, pem, is defined in Eq. (fTZt . Again we 
have neglected the factors related to stimulated emission. 

For Eq. ( fTSI l we have assumed that the emission profile for 
every new electron that was added to the 2p state is given by 
Eq. ( fT3T l, regardless of whether the electron came from the con- 
tinuum or from some excited state. In the absence of collisions 
(a very good approximation for the expanding Universe) this is 
the standard approach, in which fresh electrons, i.e. those that 
have not reached the 2p state by a line scattering ev ent, lead 
to a natural excitation of the 2p state (e.g. see p. 433 Mihalas' 
11978 ). Note that this also implicitly means that any transition of 
electrons from the 2p state to higher levels effectively leads to a 
complete redistribution of photons in the Lyman a line. 

One does expect some corrections related to these approxi- 
mations, since even for the real line emission process the history 
of the electron should mat ter (e.g. due to two-photon processes 
IChluba & Sunvaevl2008bl) . However, this problem is beyond the 
scope of this paper. 

2.2.4. Line absorption term 

In the standard formulation (e.g. see iMihalasll 19781 p. 278) the 
profile for real line absorption is usually assumed to have the 
same shape as the natural emission profile (T3[ . In this case, us- 
ing the death probability p^, the term for real line absorption 
reads 



1 dNy 
c dt 



abs 



Ly-Q- 



: Pdhv2i Bi2Nu 



4n Avd 



Ny. 



(17) 



Here A^is is the number density of hydrogen atoms in the ground 
state, and B12 is the Einstein coefficient for Lyman a absorption. 
The idea behind Eq. (fTTl i is that only a fraction of the photons 
interacting with the Lyman a line really undergo transitions to 
higher levels or the continuum, while most of the interactions 
(;jsc = Pttm = 1 - Pi) actually should be considered as resonance 
scattering events. 

More rigorously, using the principle of detailed bal- 
ance, instead of the standard absorption coefficient o-^' - 
hv2i B12 Nu (e.g. see Mihalas 1978, p. 78), from Eq. (QSll, 
also including the effect of stimulated emission, one would de- 
duce Oy - e''[''"''2il/*^^r X af. Although especially the expo- 
nential term should lead to significant differences in the distant 
wings of the Lyman a line, we follow the standard approxima- 
tion and set ^ gi^lv-v2\]/i'Ty ^ i this expression. 

It is clear that in the distant wings other corrections also 
will become very importa nt (e.g. due to two-photon emission 
IChluba & Sunvaevll2008bh . but a full consideration of these as- 
pects is beyond the scope of this paper. However, in the stan- 
dard formulation, i.e. setting ^ e'^^^^^^'^^''^^ ~ 1, akeady at 
I-^dI k, 100 - 1000 a blackbody distribution is not exactly con- 
served in full equilibrium. At the level of accuracy required in 
the cosmological recombination problem this aspect will have to 
be resolved. 



2.2.5. Final line emission and absorption term 

With Eq. ( fTSI l and ([TtI i one can now write down the collision 
term for real line emission and absorption as 



C[N 



vJle/ 



C/>(V) 



a = Pern ^2p ' PA CTi Nu (f) Ny 



= PdO-rNu(f>(v) {N^m-N 

R 



N,, 



2p 



2^21 gis 



R 



2p 



Pd hV2iBi2Nu c2 g2p R2^Nu 

Here we used the resonant scattering cross section 

^ _ hv2i B\2 ^ 7Te^ fi2 ^ 
47r Avd nigC Avd 

(l+z) 



(18a) 
(18b) 



1.91 X lO^'^cm^ 



JiA^^a 
2 

-1/2 



(19) 



Ihvl, 



and the Einstein relations A21 = —4^ B21 and B21 = fr^i2, 
where /I21 = c/v2i is the Lyman a wavelength, and /12 is 
the absorption oscillator strength of the Lyman a transition. 
Note that A^em is only a function of redshift, but not frequency. 
Furthermore, due to the factor R^^ it depends on the solution for 
the population of the higher levels. 



2.3. Transfer equation including line emission, line 

absorption and coherent scattering In the lab frame 

For coherent scattering in the lab frame no redistribution of pho- 
tons over frequency occurs. Using Eq. O and ( fTSl l. the time- 
dependent transfer equation therefore reads 



1 dN^ 

c ~dr 



Pd o", 



,Nu<t>{v)[N,^-N, 



(20) 



with V - x{ \ + z), and .^Vem = N^m/i 1 +zf', where A^em is defined 
by Eq. ( I18bl ). For the initial condition Nx(Zs) - Nf, where Zs 
is a redshift well before the epoch of hydrogen cosmological 
recombination, this equation formally has the simple solution 



NAz) = Nf - 



[N.n>iz)-m']d,'e 



-Tabs(-V,z',Z) 



dz'. (21) 



Here T^t,^{x, z',z) is defined by 

co-^Nu 



-f 

Jv 



PA 



PA 



HH+z) 
co-.Nu 



H 



(x[l+z])dz 



(v) — . 

V 



(22a) 



(22b) 



In Eq. ( I22bl ) we have used the substitution v = x(l + z), so that 
the cuiTent redshift can be found from 1 -Hi - v(l + z)/v. Note 
that in the given set of variables Nf does not explicitly depend 
on redshift, so we omitted it in the notation. 

Returning to physical coordinates one can finally write 



AN:°\z) = INUz) - Nt (z)]j®"'\z') d.e 



0™n(/) = _ 



^em(z) - Nt, 



where X2\ - V21 /(I +z) and in expression 
use X - vl{\ + z). 



-T*,(v,z',z)jj^/^ (23a) 
(23b) 

for Tabs one should 
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2.4. Transfer equation including line emission, line 
absorption and complete redistribution 

In the case of complete redistribution one has to add the term 
(see e.g. .Mih alas 1978) 

C[A^v]lr ^ P.C (r,Nu<f>(v) [n - Ny] (24) 

to Eq. dlOll. Here = / ^A^^ dvdQ. This then yields 



1 dN, 

c dt 



(25) 



with Nf^ = PdNem + Psc^, and N = N/{1 + z)^. Although 
complete redistribution is not appropriate for conditions valid 
in the expanding Universe (practically no collisions), it is used 
many times in the literature, in particular for the derivation of 
the Sobolev escape probability. One should mention that in this 
approach it is assumed that even those photons scattering in the 
very distant red wing of the Lyman a resonance can directly re- 
turn to the line center in one scattering event. With Doppler re- 
distribution, wh ich is described by so-call type-II redistribution 
(lHummerlll962h . this is only possible after many scatterings (if 
at all), or by a real line absorption event. 

Comparing Eq. ( |25l ) with Eq. ( |20] l. in physical coordinates 
one can directly write down the solution as 



wcr ( >\ _ ATI 



(26a) 
(26b) 



where t„ = rabslpj=i- It is clear that » Tabs since pd <K 1 
during cosmological recombination (see Fig. [Hi. 



3. The Lyman a escape problem and results for the 
escape probabilities 

In order to solve the cosmological recombination problem, the 
usual way is to separate the evolution of the photon field from 
the evolution of the matter, in particular the populations of the 
different energy states inside the hydrogen atom. This is nor- 
mally achieved using the Sobolev approximation for the opti- 
cally thick Lyman series in order to define the mean intensity 
of photons supporting the «p-state at a given time, and leads to 
the definition of the Sobolev escape probability. In this Section 
we explain the details of this approximation and compare it with 
other cases that can be solved analytically. 

3.1. The Lyman a net rate 

The net change of the number density of electrons in the 2p level 
via the Lyman a channel is given by 



dN2p 



dt 



= A21 — Nun-A2i(l+n)N2j, 

Ly-a Sis 



= A2i(l -Hn)A^2p 



gis N2p l+n 



- 1 



(27) 



where A21 = 6.27 x 10** s ' is the spontaneous 2p-ls transi- 
tion rate, A^, denotes the number density of electrons in level /. 



Furthermore we made use of the Einstein relations, and defined 
7^ with 

2hvi, 



n 



Jo 47r 



ly dvdD.. 



(28) 



According to the textbook derivation^ J k J(v2\) ~ hv2\N ~ 

hv2i j ^Ny dvdO. « with very high accuracy, since (p(v) 

is so sharply peaked at v ~ V21. 

Defining the line occupation number 



1 



- 1 



«LS Nig 

Eq. dZTb can be cast into the form 



g2p Nis 



(29) 



dN2p 



dt 



hy-a 



:A2,^A^isAnL, 



(30) 



where we have introduced An^ - n - ni^, for which we will now 
discuss different approximations below. 

3.2. Escape probability within the Sobolev approximation 

The aim is now to determine the solution for the mean occupa- 
tion number in the Lyman a resonance using the Sobolev ap- 
proximation. The two key assumptions for its derivation are (i) 
quasi-stationary evolution of the photon field and (ii) that every 
resonance scattering leads to a complete redistribution of pho- 
tons over the whole Lyman a line profile. With these assump- 
tions we can obtain the solution for the spectral distortion at 
redshift z using the results of Sect. 12.41 Under quasi-stationary 
conditions one can simply seQ &"(z') = 1 in Eq. ( |26] l, and for 
the absorption optical depth, Tcr, one has 



Tcr ( V, Z , Z) ~ Ts (Z) (p(v) dv 



r 



(31) 



where (p{v) - (p{v,z) is given by Eq. ( fT3] ), and v' = ^Tk^- 
Furthermore we introduced the Sobolev optical depth of the 
Lyman a line 



C<Ti-Nu Avd g2p A21/I21 



H 



V2l 



Nu 



(32) 



with wavelength A21 - c/v2i. For Eq. ( [3T1 ) we have assumed 
that Ts does not change significantly between z' and z. Also we 
have neglected the variation of 1/y in comparison with (f(v), and 
set 1/v X! I/V21. This approximation is normally applied in the 
literature and computation of the recombination history. 
From Eq. ( l26l ) with ©"^"^ = 1 one then obtains 

AN^-'^'iz) - [KL - Kl ] [1 - e-^'e^'"'] , (33) 
with^(v) = £ (f(v')dv'. 

* Also in the derivation of the Einstein relations this approximation 
is normally apphed. 

' In fact this approximation not only implies quasi-stationarity, i.e. 
Nl^^iz') = NSn(z), but also that one can use f/f « N^^ . 
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Figure 2. Comparison of the Sobolev optical depth, ts, as de- 
fined by Eq. ( |32] |. and the total absorption optical depth, Td - 
PdTs, for the 2- and 3 shell case. To compute ts we used the 
Recfast solution for A^is. 



3 shell Hydrogen atom, z= 11 00 



0.4 - 



0.2 




no line scattering 
complete redistribution 



10" 10 10 10 10' 10" 10 10" 10 10 



Figures. Spectral behavior of the solutions in the quasi- 
stationary approximation at redshift z - 1100. We normal- 
ized the distortion to unity at the Lyman a frequency, i.e. F = 
ANviz)/ ANv2,iz)- The solid line shows the result obtained in the 
no line scattering approximation, while the dashed line repre- 
sents the solution for the complete redistribution case. We also 
indicated the position of the Lyman /? resonance at that time. 



3.2.1. Spectral characteristics of tine solution 

As can be seen in Fig.|2] during hydrogen recombination Ts » 1. 
According to Eg. ( |33T l the photon distribution therefore varies 
strongly close tcQTs[l - x] ~ lii2, while it is basically identical 
to unityOat xd < 0. Using the wing-expansion (IA.2b of the Voigt 



* There the value of AAf"'''''(z) has decreased by a factor of 2 com- 
pared to the line center. 
' With relative accuracy better than ~ e"'"^'^. 



profile one therefore finds that this happens at a distance of about 
-u^-^. (34) 



At z ~ 1 100 one has a ~ 8.6 x 10^"^ and ts ~ 5.8 x 10^ such that 
xu ~ 2.3 X 10'''. This shows that in the complete redistribution 
approximation the photon distribution is in full equilibrium with 
the value at the line center up to extremely large distances on the 
blue side of the line center (see Fig.|3]l. 

Physically this type of redistribution does not describe the 
problem very accurately, and a much more realistic solution is 
obtained using the case of coherent scattering in the lab frame 
(see Sect. 13.31) . For example, if we consider the position of the 
Lyman /3 line at redshift z = 1 100, then in Doppler units of the 

Lyman a line one finds x^'' ~ 8000. The variation of the pho- 
ton distribution, which is important for the value of the escape 
probability (see below), occurs far beyond this value. In fact, 
xu ~ 2.3 X 10^ corresponds to about 2 times the Lyman a fre- 
quency, or 1 .5 times the ionization energy of the hydrogen atom. 
The Sobolev optical depth, ts, for conditions during recombina- 
tion is simply so large that the approximation of complete redis- 
tribution becomes unphysical. 

Furthermore, at such large distances it is even questionable 
as to why one should be able to neglect variations of the black- 
body distribution, or the factor of 1/v in the definition of Tabs- 
However, such an approximation is necessary to obtain the ex- 
pression for the Sobolev escape probability. Obviously other 
coiTections (e.g. related to two-photon processes, or the imbal- 
ance in the emission and absorption coefficient as mentioned in 
Sect. 12.2.4b will become important and even necessary to cor- 
rect for these physical discrepancies. However, as we will see 
below, in spite of all these problems the Sobolev approximation 
at the level of ~ 10% provides the coiTect answer for the escape 
probability, a fact that is very surprising. 



3.2.2. IVlean occupation number in the Lyman a and the 
Sobolev escape probability 

To obtain the mean photon occupation number in the Lyman a 
line we multiply (|33T l by tp and integrate over v. This then yields 

An = n - nP' I^p^^ n + Hem - Wvl,} [1 - -Ps] (35a) 

Pg = 1 - I Hi/^>-^s[i-;\r] 



= 1- r Axe- 
Jo 



L^— i ±, (35b) 

TS TS 



where we again have neglected the variation of nP' over the line 
profile. Here Ps denotes the standard Sobolev escape probability 
(see Fig. [1] for its redshift dependence). After some rearrange- 
ment and with n^'j w fiP^, Eq. ( I35al ) can be cast into the final 



form 



_pl ^ Pd(l - f s) , _pi, 

"P + — [Wem - nP ] 

Pi + PscPs 

Ps 



PA + PscPs 



["em - «'''] 



I 



- — [«em -«"']. 
Pd 



(36a) 
(36b) 

(36c) 



A solution similar to Eq. ([36b was also given and discussed in 
iHummer & Rvbickil (|1992|) . 
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3.2.3. Relation to the expression wliicli is normally used in 
multi-level recombination codes 



But how does Eq. ( |36] | actually relate to the expression 



nP'] = nP' + [nL-nP'](l-Ps) (37) 



that is normally used (cf. lSeager et al.ll2000h in computations of 
the hydrogen recombination problem? To understand this con- 
nection the key ingredient is the quasi-stationary solution for the 
2p population. In fact this approximation should always be pos- 
sible, even if the spectral evolution is non-stationary, simply be- 
cause the re-adjustment of the 2p population after some changes 
in the spectrum is so fast. 

With Eqs. ( fTOb and ( l30b , the rate equation governing the time 
evolution of the 2p state can be cast in the form 



dX 



2p 



dt 



2lA^2p 



- 1 



(38) 



3 shell Hydrogen atom 




''600 



800 



1000 



1200 



1400 



1600 



1800 



2000 



where we directly neglected induced terms, and introduced 
X2p = N2p/Nu, where A^h denotes the total number density of 



hydrogen nuclei. With 



one readily finds 



Figure 4. Relative difference between the escape probability 
'''' and the standard Sobolev escape probability, Ps- 



l-Pcmi'd 

The death probability for a 3 shell hydrogen atom was used 



A21+«2p-A2,iL 



Pd»L _ 
"L - Pemii R 



Inserting this into Eq. ( llSbl i therefore yields 



(39) 



(40) 



If we now use this in Eq. ([35j one immediately finds «ci = "st- 
Therefore the result (|36] l is completely equivalent to Eq. ( [37T i. 

However, in one case factors are expressed in terms of nem, 
while in the other case «l is used. From Eq. (I36bb and (iJTl i with 
"ci = «st one can easily show 



5Pl 



1 + '- [nL-nP']- 

Pd 



(41) 



This implies that nem ~ "l- However, since in the Lyman a rate, 
Eq. ( l30t , the main term (~ «l) cancels, the small difference «em- 
-[hl - nP'] cannot be neglected. Replacing the first rism 



in Eq. ( |36ct with this expression and neglecting terms C)(f s / Pd)^ 
one can directly recover Eq. dJTl i. However, here it is important 
that in Eq. (|36c| i there is partial cancellation of terms 0(Ps/pd) 
from the first and second «em, leaving a much smaller residual 

~PS[«L-«P']. 



3.3. Escape probability for tine case of coinerent scattering 

In the absence of line scattering, or equivalently for coherent 
scattering in the lab frame, the solution of the transfer equation 
is given by Eq. (|23] |. Under quasi- stationary conditions (and with 
A^P' X A^P'j) one again has &^°^{z') - 1, and also it is possible to 
use T^^j, ^ Pd Tcr, where is defined by Eq. OTt . Then one can 
write 



(42) 



with Td = PdTs. 



3.3.1 . Spectral characteristics of the solution 

Looking at Fig. |2]it is clear that Td <K ts at all relevant redshifts, 
so that hNl"^''^^{7) should change strongly much closer to the line 
center than in the complete redistribution case, Eq. ( |33] |. If we 
again want to estimate where the photon distribution ( l42b varies 
most rapidly, assuming that this happens in the blue wing of the 
Lyman a line, we can find 



fl PdTs 

;7r ln2 



(43) 



At z ~ 1 100 one has pd ~ 10^'^, such that with Eq. ( [34l l one has 
xd ~ 30 (see Fig. [3). This result is much closer to the solution 
that would be obtained when using the more realistic type-II re- 
distribution for the Lyman a resonance scattering process, where 
in the quasi-stationary approximation the photon distribution for 
typical conditions in our Universe strongly varies at distances of 
a few hundred Doppler width (e.g. see Rvbic ki & deirAntoniol 
,1994). We already checked this point and found very similar re- 
sults dChluba & Sunvaevll2009l) . 



3.3.2. Mean occupation number in the Lyman a and the 
escape probability 

With ( l42l i and the same simplifications that were mentioned 
above in connection with Eq. ( l36b one then obtains 



Wcoh 



1 - e-^" T' 



Td 



nP'] 

Pd' 



PdTS 



(44a) 
(44b) 



Note that Pd is very similar to the standard Sobolev escape prob- 
abiUty, Ps - [1 - e^'^^]/Ts, with the only difference that in 
general Td < ts and hence Pd > Ps (cf. Fig. [T] and |2]i. Also 
one can directly see that in lowest order Eq. ( l44l i is identical to 
Eq. ( |36c| i. This already suggests that in both the complete redis- 
tribution and the no line scattering approximation the answer for 
h is nearly the same, with differences of the order 0{Ps / Pd)^- 
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Again using the quasi-stationary solution for the 2p popula- 
tion, we can replace nem applying the expression Eq. ( |40l ). Then 
solving for n one finds 

PdPd r -Pin fAC\ 

«coh ~ "L - -; ;r["L - n ]. (45) 

1 - f em ^^d 

Comparing this with the standard form Eq. (|37] |. it is again clear 
that for PemPd ^ 1 and rs s> Td s> 1 one has hcob ~ «st- 
Looking at Fig. [1] shows that this should be the case at most of 
the redshifts relevant for the recombination of hydrogen, and that 
the differences between the complete redistribution and no scat- 
tering case should not exceed the level of ~ 10"^ in the redshift 
range 800 < z < 1600. 

In Fig. ID we present a more detailed comparison and in- 
deed find practically no important difference to the standard 
Sobolev case. This result is somehow surprising, since the as- 
sumption of complete redistribution leads to a totally different 
(and physically unrealistic) solution for the photon distribution. 
Still, the final result is comparable. This is due to the fact that 
the changes in the shape of the photon distribution are compen- 
sated by changes in the amplitude of the spectrum close to the 
line center, as already explained in connection with Eq. ( |4TI ). 

Note that in the case of 2 shells the differences would be 
much greater, since one can find f d ~ 1 at redshifts relevant 
for recombination (see Fig.[T]i. This again shows that one has to 
include at least 3 shells in the computation, in order to obtain 
meaningful results. 

3.3.3. Escape probability in tine limit pa 1 

It is also illustrative to look at the solution in the limit pd ^ 1 . 
Physically, in the current formulation of the problem this should 
give the same answer as in the approximation of complete redis- 
tribution. This is because for pi - \ every electron entering the 
2p level via the Lyman a channel will pass through the contin- 
uum or some higher shell, where it will forget its history. It will 
be replaced by another fresh electron, with a natural line pro- 
file, as in the complete redistribution approximation for a line 
scattering event. 

From Eq. ( |45] |. with pd = 1, Psc - and = Ps, it is quite 
obvious that Wcoh = n„, but can one also see this directly from 
Eq. ( |44a| ), which in the first place only leads to Wcoh - "em - 
^'s[nem - «''']■ Here apparently «em = "l, a result that indeed can 
be confirmed with Eq. ( |40l ), so that also «coh = «ci follows. 

3.3.4. Until what distance from the line center is the shape of 
the photon distribution important? 

The escape probability, P^, was obtained from the integral over 
the Lyman a line profile. If we only integrate up to some fre- 
quency Vm, then one has 

^'d(Vm) = (fiiv) e^^'-e^''^^"' dv = 

Jo Td 

g-Td[l-A-m] 

* , (46) 

Td 

with Xm - xi'^m)- Therefore the relative difference in the value 
of Pd is given by 

APd(v,n) Pd(v,n) - Pd 



If we assume that Vm > V2i is already in the damping wing then 
with the approximation Eq. jA.H of the Voigt profile one obtains 



APd(Vm) 
Pd 



a Td 

TTXu' 



(48) 



-16% [^] ', so that in the 



At z = 1100 this yields 

no-scattering approximation for ~ 10%, ~ 1%, and ~ 0.1% ac- 
curacy one has to know the spectrum up to xd ~ 10^, jcd ~ 10^, 
and Xd ~ lO'*. This shows how important the knowledge of the 
solution for the photon distribution in the distant wings is. 

In the case of complete redistribution one can easily show 

that = -16% [j^Y^ at z = 1 100. This implies that for ~ 

10%, ~ 1%, and ~ 0.1% accuracy one has to know the spectrum 
up to Xd ~ 10^, Xd ~ 10^, and xd ~ 10^. Let us emphasize again 
that these are extremely large (even unphysical) distances from 
the Lyman a resonance. However, it is in these regions where 
the value of the Sobolev escape probability is formed. 



3.4. Effective escape probability using ttie time-dependent 
solution 

With the solution (l23T l we can also describe the time-dependence 
of h within the approximation of coherent scattering in the lab 
frame. Although one does expect some modifications when ac- 
counting for partial freq uency redistribution, our computations 
(IChluba & SunvaevI l2009i) show that the additional correction 
will be dominated by the influence of line reco iP^ which has 
been addressed in iGrachev & DubrovichI (1200^ . However, the 
time-dependent correction that is considered here turns out to be 
much larger, so that we shall focus on this only. Below we now 
provide a detailed discussion of the time-dependent correction 
in the case of coherent scattering in the lab frame, introducing 
an effective escape probability, which then can be used in com- 
putations of the cosmological recombination history. 



3.4.1. Escape probability during recombination without 
redistribution but with full time-dependence 

Using the time-dependent solution for the case of no redistribu- 
tion, Eq. ( l23T l. it is possible to write 

2v2i Jo 

= nP' + An,M) J ^(v) dvj &'"\z') d^e-^^^^--^ dz' 
= nP' + A«e™a)(l-0 (49a) 
Pl^^l+ J^^ifiiv) dv J" 0'=°''(z') S.-e-^^^"'^'-^^ dz', (49b) 



where we have Anem(z) = We 



Here it is very important to mention that one has to use 



0"^° (z) as defined by Eq. (I23bt but evaluate the blackbody dis- 
tribution at the line center only, i.e. use &'^°^{z') ~ [Nemiz') - 
]/[Nem(z)-Nf^,] with xl. = V21 /[I H-z']- This is necessary in 
order to be consistent with the formulation of line emission and 



rdfl-^ml«l 



-rd[l-^..]_l 4 (47) 



Including atomic recoil we find a correction of AP/P ~ 4% at 
z ~ 1100 and AP/P ~ 6% at z ~ 800 t o the Sobolev escape probabil - 
ity, which, in reasonable agreement with lGrachev & DubrovichI i2008l) . 
leads to ANJN^ ~ -1.2% at z ~ 950. 
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absorption processes, which, as mentioned in Sect. 12.2.41 in full 
equilibrium does not exactly conserve a blackbody distribution. 
A more consistent formulation will be given in a future paper, 
but the result for the pure time-dependent correction should be 
very similar 

Eq. ( 149 a| ) provides the time-dependent solution for h(z), 
when the ionization history is known until z- However, in real 
calculations Eq. ( |49a| i is not very useful, since the evaluation of 
the integral is rather time-consuming. With Eq. ( |49b| i we defined 
an effective escape probability, which can be compared with the 
result in the full quasi-stationary case. The differences will be 
due to non-stationary contributions in the evolution of the pho- 
ton distribution, and can be iteratively used in computations of 
the recombination history. Since the correction is expected to be 
small, even the first iteration should give a rather good answer. 

To obtain the difference from the Sobolev escape probability, 
one again has to use the quasi-stationary solution for the 2p state, 
leading to relation (l40t . With this one can eliminate Wem from 
( |49] l, and bring the expression for n in the standard form ( |37] |. 
This yields 



3 shell Hydrogen atom 



r^coh 



(z) = nL-/^esc ["L-n"'] 



esc 



1 - P.cP[, 



(50a) 
(50b) 



Now can be directly compared with the Sobolev escape 
probability. 

Looking at ( |49bt it is clear that there are two sources for 
the time-dependent correction. The first comes from the time- 
dependence of 0'^°*', while the second is due to modifications in 
the absorption optical depth. Tabs- Below we now discuss each 
correction separately. 

3.4.2. Neglecting the time-dependence of 0™'' 
If we set 0'^°'' = 1 in Eq. ( |49b| ) one obtains 

pZ(z)^ ^(v)e-^-('''--^)dy. (51) 
Jo 

With this expression it is possible to take into account the time- 
dependent corrections that are only due to the modifications of 
Tabs in comparison to the quasi-stationary case (see Sect. |3.3l l. 

First, it is clear that due to the v-dependence of the absorption 
cross section the total absorption optical depth depends strongly 
on the initial frequency of the emitted photon. For example, if 
a photon is emitted on the blue side of the Lyman a resonance, 
then after some redshifting it will come close to the Doppler core 
of the Lyman a line, where it will be absorbed with extremely 
high probability. Depending on the initial distance to the Doppler 
core, this will take some time, during which the properties of the 
medium (e.g. the ionization degree) may have changed signifi- 
cantly. Similarly, photons emitted in the very distant red wing of 
the Lyman a line may redshift for a very long time, before they 
will be reabsorbed, if at all. 

At high z the total absorption optical depth is expected to 
mainly vary due to the changes in the number density of ionized 
hydrogen atoms, and at low redshifts because of the steep drop 
in Pi. If for given initial frequency v' of an emitted photon the 
time it takes until this photon is reabsorbed (Tabs ~ 1) is similar 
to the Hubble time, then these changes may be important. 

If the considered photon was emitted close to the Lyman a 
line center, the absorption optical depth is dominated by its value 
inside the Doppler core, where photons only travel a very short 
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Figures. Differences between the escape probability P = 



1-p.cPi; 



where P^^iz) is given by Eq. ( BTT i. and the standard 

Sobolev escape probability, Ps. We used the death probability 
for the 3 shell hydrogen atom. 



distance (a small fraction of the Doppler width), before being 
reabsorbed. In this case, the quasi-stationary approximation cer- 
tainly is valid with very high accuracy, since Av/v <K 1 between 
emission and absorption implies Az/z <c 1, so that the medium 
has not changed very much. However, when the photon is ini- 
tially released in the distant red or blue wing of the Lyman a 
resonance, it can redshift for a much longer time before being 
reabsorbed, so that changes in the medium, in particular the ion- 
ization degree and death probability, may play an important role. 

In Fig. |5] we show the direct comparison of the escape prob- 
ability that follows from Eq. ( fSTT i with the Sobolev escape prob- 
ability. At very low and very high redshifts the correction due to 
the pure time-dependence of Tabs becomes very small. The dif- 
ference that is seen close to z ~ 600 and z ~ 1800 is only related 
to the correction coming from the coherent scattering approxi- 
mation (see Fig.HJi. In both cases this behavior can be explained 
by the fact that the importance of the wings for the total value of 
the escape probability decreases. Photons escape directly from 
the Doppler core, so that the contributions to the value of PeVn(z) 
can be considered quasi-stationary. 

To understand the behavior at intermediate redshift, it is im- 
portant that before the maximum of Td around Zmax ~ 1300 (cf. 
Fig. |2]i, one expects that independent of the considered fre- 
quency, Tabs(v, z',z) is smaller than in the quasi-stationary ap- 
proximation t'^Is(v, z' , z) - Td(z) <p(v")dv". This is simply 
because at z' > z > Zmax the value of Td(z') < Td(z). At those 
times the time-dependent modifications of Tabs should therefore 
result in a positive correction to the effective escape probabil- 
ity (cf. Fig. Is). With a similar argument, at redshift z < 1300 
the correction in the escape probabihty should be negative, as is 
seen in Fig.|5] 



" This statement is not completely correct, since in Eq. l l22bt we do 
take into account the factor 1/v. However, this only affects the very 
distant blue wing. 
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Figure 6. The spectral distortion AA^v far away from the Lyman a line center at different stages of hydrogen recombination. In all 
cases redistribution of photons over frequency was neglected. We normalized the distortion to unity at the Lyman a frequency, i.e. 
F = ANy{z)/ ANv^_j{z)- The dotted curves give the spectral distortion in the quasi-stationary approximation, F''^(v) = 1 - e^'^'^e'^'^^'^K 
For the dashed curves we neglected the variation of 0'^°'' in the solution ( |23] |. leading to Fo(v) - I - e^'^'->"^^'~''~\ with Zs — » The 
solid lines show the spectral distortion in the time-dependent case, i.e. F{v) = Fo{v) + AF(v), where AF{v) is defined by Eq. ( |52bt . 
For comparison, the dash-dotted curves give the spectral distortion in the Sobolev approximation, F^(v) - 1 - e^'^^ e^^^'^^K We also 
indicated the position of the maximum of the redshifted Lyman a distortion, which appeared at redshift z ~ 1400. 



3.4.3. Correction due to the time-dependence of 0' 



,coh 



In Sect. 13.4.^ we have neglected the time-dependence of 0'^°*'. 
This factor describes how much the photon emission process 
varies as a function of time, which in the present approximation 
is independent of frequency (see comment in Sect. |3.4.TT i. 
With Eq. ( fSTI l and ( |49bt one can define 



APl^ = ^ ^(v) dv £ [&'°\z') - l] a-e-"'*^*"-' ^> dz' 

= - [ tp(v)AF(v)dv 
Jo 

AF(v) = £ [®'°\z') - l] 5,,e-^*<''-"'-"'dz'. 



(52a) 
(52b) 



With this expression it is now possible to calculat43 the addi- 
tional correction to the escape probability coming from the vari- 
ation of 0'^°'^ over time. 

In order to understand the final result we first con- 
sider the behavior of the inner integrand (15 2 al l at different 
stages of hydrogen recombination. Since the function F(v) = 
j^' 0coh(-2') ()_,g-TL,b,(v,z',z) (j^' is identical to the Lyman a spectral 
distortion in the no redistribution approximation, but normalized 
to its value at the line center, i.e. F{v) = ANv(z) / ANv2i(z), it is 



The evaluation of the integral l l52at is rather cumbersome. It is most 
important that for a fixed frequency v at redshift z the inner integral 
varies most strongly at z' ~ max(z, V2i(l -I- z)/v). 
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Figure 7. The spectral distortion ANy close to the Lyman a line center at different stages of hydrogen recombination. In all cases 
redistribution of photons over frequency was neglected. We normalized the distortion to unity at the Lyman a frequency, i.e. F = 
ANviz)/ ANv,t(z) and show the difference to the quasi-stationary solution, F'^^{v) = 1 -e"'"''e'"'J^*'''. For the dashed curves we neglected 
the variation of 0'^°'' in the computation of the escape probability, while for the solid lines the spectral distortion in the time- 
dependent case was taken into account. 



illustrative to define 

f^'iv) = 1 - e-^^e^i*W (53a) 

Fo(v) = 1 - g-^^i.^t"'^-^) (53b) 

F{v) = Foiv) + AF(v), (53c) 

in addition to Eq. (|52b| i. Here F'^^{v) represents the Lyman a 
spectral distortion in the full quasi-stationary approximation, for 
F()(v) only the variation of 0'^°*' is neglected, and F{v) is the 
time-dependent distortion in the no redistribution approxima- 
tion. From these functions the corresponding escape probabil- 
ities can be obtained by Pgsc - ^ ~ J dv, which then 
can be inserted in Eq. (|50b| i when comparing with Ps. Therefore 
it is clear that if F{v) < Fi^(v) at all frequencies, the resulting 
effective escape probability should be slightly larger than Pj- 

In Fig.|6]and|7]we illustrate the behavior of the functions ( l53T l 
at different stages of hydrogen recombination. We used the solu- 



tion for the populations in the 3 shell case as given by our rn ulti- 
level code jRubino-Martm et al.l 120061: IChluba et al]|2007h . As 
expected, in all cases F''^(v) and Fo(v) are very close to unity at 
jcd ^ and then drops very fast toward zero at xd > 0. Also 
Fig. |2] clearly shows that Fi^(v) ^ Fo(v) in the red wing and 
the Doppler core of the Lyman a resonance. This is expected, 
since at v < V2\ always Tabs » 1, so that its exact value does 
not matter. Furthermore we can observe a change in the sign of 
the difference Fq{v) - F''^(v) in the blue wing when going from 
high to lower redshift. At z = 1400 one can clearly see that 
Fo(v) < Fi'(v) in the range < xd < 100, so that Pe° (z) > Pa 
is expected, in agreement with the results presented in Fig. |5] 
On the other hand, in all the other cases shown F(){v) > F''^(v) 
at < jcd < 100 so that one should find Pem(z) < Pd, again 
confirming the results given in Fig.|5] 

Because of the steep drop of F'*^(v) and F()(v) at a few 
Doppler width above the line center, the main contribution to 
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the escape probability clearly comes from rather close to the line 
center. However, at the level of percent the shape of the distortion 
up to a few hundred or thousand Doppler widths is important. 

If we now look at the spectral distortion in the time- 
dependent approximation, F{v) - Fo(v) + AF{v), we can see 
that at all stages the variations of 0'^°*' with redshift become im- 
portant outside the Doppler core. From Fig.|6]we can distinguish 
in more detail the following regimes: (i) at redshifts z > 1400 the 
distant wing distortion is smaller than in the quasi-stationary ap- 
proximation. This is because at redshifts much before the time 
under consideration the emission in the Lyman a transition was 
very inefficient, so that until then not many photons can have ap- 
peared or reached large distances from the Lyman a line. The 
slope of the red wing distortion is positive close to the line cen- 
ter; (ii) at redshifts z < 1400 the distortion in the blue wing and 
nearby red wing is greater than in the quasi-stationary approx- 
imation. The production rate of Lyman a photons has already 
passed its maximum (at z ~ 1400), so that at the current line cen- 
ter/ewer photons are produced the lower the redshift becomes. 
The slope of the red wing distortion is negative close to the line 
center 

It is also clear that in case (i) the value of An is smaller than 
in the quasi-stationary approximation, while it is expected to be 
larger in case (ii). According to the definition ( |49] l this implies 
that in the former case the effective escape probability is higher 
than in the quasi-stationary approximation, while it is lower in 
the latter case. In Fig. [8] we can see that these expectations are 
true (see solid line). The total correction due to excess or a deficit 
of photons leads to a total decrease of the effective escape prob- 

abihty at z < 1400 that reaches AP/P 5.8% at z ~ 1140, 

while it results in an increase of AP/P — 1-3.7% atz~ 1510. 

Although in the escape integral the distortion in the vicin- 
ity of the line center mainly contributes, at the percent level the 
distant wings are also important. As we have seen in Fig. |6] the 
red wing distortion due to the Lyman a transition can exceed 
the distortion close to the line center by a large amount. In this 
case the question is how much the very distant wings actually 
contribute to the total correction shown in Fig. |8] For this we 
computed AP[.^(z), but excluding the correction at xu < -10^. 
Looking at the boxed curve in Fig. |8] shows that the very dis- 
tant red wings contribute about AP/P ~ 0.6% at z ~ 1500, and 

AP/P 1.3% at z ~ 1100. This is an important point, since 

in the very distant wings other processes related to the formu- 
lation of the problem will also become important (i.e. due to 
changes in the absorption profiles, when considering the prob- 
lem as a two-photon process), so that one expects additional re- 
visions for contributions from the very distant wings. However, 
here the corrections mainly seems to come from regions in the 
vicinity (IjcqI < 10^ - lO^*) of the Lyman a line center. 

At low redshifts (z < 800 - 900) one can observe an addi- 
tional strong decrease in the effective escape probability. This is 
due to the additional re-excitation of electrons by the distortion 
on the blue side of the Lyman a resonance. Most of the photons 
in this part of the spectrum have been emitted much earlier, at 
times around the maximum of the Lyman a emission (z ~ 1400). 
This also explains the huge difference to the quasi-stationary so- 
lution: as one can see in Fig.|6] at z < 800 the amount of photons 
exceeds the spectral distortion obtained in the quasi-stationary 
approximation by about two orders of magnitude. The spectral 
distortion is only a factor of ~ 100 below the emission in the line 
center Looking at Fig. [T] very close to the Lyman a line center 
some differences also are visible, which at the percent level do 
matter. 
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Figures. Correction to the Sobolev escape probability due to 
variations of 0'^°'^. For the solid curve we used AP^^iz) accord- 
ing to Eq. ( i52l ). The boxed curve was computed with Eq. ( i52l i. 
but including the variation of 0'^°'^ only at xd > -10^. The other 
two curves were obtained setting 0'^°*' = 1 in the red wing (stars) 
or blue wing (crosses), respectively. 

To show that the distortions on the blue side of the Lyman 
a are responsible for this re-excitation we also computed the 
correction only including the non-stationary contributions for 
the red side, but setting 0'^°'^ = 1 for evaluations on the blue 
side. The result is also shown in Fig. |6] (stars). For complete- 
ness we also gave the curve when only including the corrections 
on the blue side of the line. As one can see, at z > 1100 the 
red and blue wing corrections are very similar. However, at low 
redshifts the blue wing correction clearly dominates, supporting 
the statement made above. Again one can expect some changes 
in the conclusions when treating the problem in the full two- 
photon formulation, since the emission of photons at high fre- 
quencies will be significantly less in the two-photon treatment, 
simply due to the fact that due to energy conservation the emis- 
sion pro files do not extend to arbitrarily high frequencies (e.g. 
see Chlu ba & Sunvaevll2008bl) . However, the corrections to the 
escape probability at z < 800 do not propagate very strongly to 
the ionization history, so that here we do not consider this un- 
physical aspect of the solution any further. 

3.4.4. Dependence of the effective escape probability on tine 
included number of hydrogen shells 

Although the death probability pi does not dependent on the so- 
lution of the recombination problem, the amount of fresh elec- 
trons injected into the Lyman a line depends on the populations 
of the excited states. Therefore the strength of the Lyman a line 
strongly depends on the total number of shells that are included 
(Rubifio-Martm et al. 2006). Here in particular the low-redshift 
tail (z < 800) will be affected, and hence there one also expects 
changes in the correction to the effective escape probability. 

In Fig. [8] we show the differences in the escape probability 
when including more shells. We used the numerical solution for 
the excited levels as obtained with our multi-level hydrogen code 
dChluba et al.ll2007l) . At redshifts z > 1200 the result is practi- 
cally unaffected by the total number of hydrogen shells that are 
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Figure 9. Total correction to the Sobolev escape probability. All 
curves were computed using the time-dependent solution ac- 
cording to Eq. ( |49bt . but including a different number of shells. 



included. In particular the result seems converged when includ- 
ing ~ 4-5 shells, leading to a total corrections of AP/P 7.6% 

at z ~ 1 100, while it results in an increase of AP/P — 1-6.4% at 
z ~ 1490. At redshifts 900 < z < 1200, some small changes 
are still visible when including more than 3 shells, but again 
the result seems to remain unchanged when including more than 
~ 4 - 5 shells. For the computation of the CMB power spectra 
the corrections in this range are most important (see Sect.HJi. 

At lower redshifts, however, the result still changes notably. 
At z ~ 600 the correction increases by about 6.4% when includ- 
ing 5 shells, and for 10 shells even by about 10%. This can be 
explained when realizing that the total emission in the Lyman a 
line at low redshifts becomes less when including more shells. 
Therefore the wing emission from redshifts around the maxi- 
mum of the Lyman a emission (z ~ 1400), which practically 
remains unchanged, becomes more important, being able to re- 
excite the 2p state as explained in the previous paragraph. As we 
mentioned akeady this aspect will probably be affected when 
including corrections to the emission and absorption profiles ac- 
cording to the two-photon formulation. Furthermore, as we will 
see in the next Sect, this low redshift tail is not so important for 
the predictions of the CMB power spectra. 



3.4.5. Dependence of the effective escape probability on tine 
distance from tine line center 

As a last point, we want to answer the question of where in the 
case of APg^ the main coiTection actually comes from. We have 
already seen in Sect. 13.4.31 that the very distant red wing con- 
tributes at the level of percent. Also we have seen that at red- 
shifts z > 800 - 900 the blue and red sides of the Lyman a line 
give similar contributions, while at low redshift due to the self- 
feedback the blue wing clearly dominates. 

In Fig.[TO]we show the total coiTection to the Sobolev escape 
probability when only including the time-dependent connection 
for a given central region around the resonance in the computa- 
tion of APgi^. It is obvious that the innermost Doppler core (+1 
Doppler width) does not contribute much to the result. This is 
expected, since there the photon distribution should evolve as 
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Figure 10. Total correction to the Sobolev escape probability 
for the 5 shell hydrogen atom. All curves were computed us- 
ing the time-dependent solution, - f [.'^ + AP^^ accord- 
ing to Eqs. ( I49bb . ( BTl i and ( |52| |. The term Pgi^ was fully in- 
cluded, but for the contribution from AP^^ the time-dependent 
correction was only taken into account for a given central region 
|xd| < const around the resonance, as labeled respectively. 



in the quasi-stationary approximation, with very high accuracy. 
This fact can also be seen in Fig.|7] where close to the line center 
the deviation of the photon distribution for the quasi-stationary 
solution is very small. 

At low redshifts (z < 800 - 900) the region up to ±4 to 
+ 10 Doppler width seems to be quite important. As explained 
in Sect. 13.4.31 there the coiTection is mainly because of self- 
feedback, which is strongest where photons really are reab- 
sorbed. However, at practically all other redshifts one can clearly 
see that the distant wings contribute significantly. Even within 
±10^ Doppler width the deviations of the spectrum from the 
quasi-stationary solution are important. 

3.5. Dependence of the escape probability on tlie sliape of 
fhe emission and absorption profile 

As mentioned in the introduction, in the Sobolev approximation 
it is well known that the result for the escape probability does not 
depend on the shape of the Lyman a emission profile. Looking 
at the derivation of expression ( I35bb for Ps it is clear that in 
addition to the condition of quasi-stationarity one needs 0em = 
<p!Lhti = Ipse, i-e. the equality of the line emission, line absorption, 
and line scattering profile. This conclusion is also reached in the 
case of no line scattering (Sect. I3.3l l leading to as given by 
Eq. ( |44bl ). 

However, if 0em ^ <f>dhs then the situation is more compli- 
cated. Starting with Eq. ( |20] i, but allowing 0em 9^ (pabs, one can 
find the solution 

AA^:(z) = [NUz) - Nt (z)]jy(z') d,e-''-^^''''''''> dz', (54a) 

(54b) 
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Figure 11. Relative change in the number density of free elec- 
trons when allowing for a constant relative change in the Sobolev 
escape probability of the Lyman a transition. At maximum 
ANJN, ~ -lAP/P. 



with V.' = V j^, and rabs(v, as defined by Eq. (l22l l but re- 
placing (f) — > 0ab,s- With this, the mean occupation number in the 
Lyman a line can be directly computed leading to 



3 shell Hydrogen atom 




v .'■ total time tleptindcnt con"cction - 

-1 — Griichev & Dubrovich — 

(2008) 



-1.2 - , - ' sum of both — 

.1 4L I I I T \ , I , \ , \ , I , \ ^ 

600 700 800 900 1000 1100 1200 1300 1400 1500 1600 

Figure 12. Estimated relative changes in the number density 
of free electrons when including various physical processes. 
The curves were obtained by simply multiplying the computed 
change in the escape probability for the 10 shell case as given in 
Fig.lSby the curve in Fig.[TI]for AP/P ^ + \% ( thick soHd fine). 
We also show the result of Grac hev & Dubrovich. (.2008.) and the 
resulting sum of both (dashed). 



(2 

^™h(2) = ^Pl + _ ^,bs(v)A<(z)dv 
2^21 Jo 

= «P' + Anem(z) Vabs(v) dv J &\z') d.^e-^-'^^'"''' ''^ dz' 
= «P' + A«em(z)(l-PL) (55a) 
^em = 1 + j| ^abs(v) dv J ®\z) c^^-e-^**''--^^ dz', (55b) 

with the same definitions as in Eq. (|49T l. Note that here one has 
to compute h using the absorption profile function. 

Now it is easy to show that for 0^ = 1, Tabs = Pd^sX^ and 
<pem = f^abs - 4>^ One directly obtains the solution Eq. (l44l l. Here 
the crucial point is that it is possible to introduce the variable 
Xiy) - ^ 0(v')dv', so that there is no direct dependence on 
(j). However, it is clear that already for 0em ^ 0ab,s, in general 
0^ ?i 1, implying that this coordinate transformation can only 
be achieved approximately, so that the result will depend on the 
shape of the Lyman a line. Similarly, in the case 0em = 0ab,s, but 
for ©" 1 and Tab.s f{z)x^ where f{z) is a function of red- 
shift only, the result for the escape probability will depend on 0. 
In that case, it will be important how large the deviations from 
quasi-stationarity are in the range, where most of the contribu- 
tions to P originate. If for example the profile is a pure Doppler 
core, then any photon that is emitted will at most travel, scat- 
ter and redshift over a characteristic length Av/v ~ few x 10"^ 
before escaping. This corresponds to Az/z ~ 10"^, so that the 
properties of the medium have not changed very much, and the 
correction should be ANs/Ns ~ few x 10"-'%. However, if, as 
in the real problem, radiative transfer is occuring in the distant 
Lorentz wings, these coiTections will be important. 



4. Corrections to the ionization history 

With the results of the previous Section it is possible to esti- 
mate the expected changes in the ionization history. In Fig. [TT] 
we show how a constant difference in the Sobolev escape proba- 
bility affects the ionization history. We can see that the response 
is roughly proportional to the given AP/P. Therefore we can use 
the curve for AP/P = 1% to estimate the changes in the ioniza- 
tion history for the results given above. Since all the corrections 
are small, one expects a small additional correction, when com- 
puting the escape probability for this modified ionization history. 
For the purpose of this paper this approximation is sufficient. 

Fig. [TT]also shows that percent level corrections to the es- 
cape probability do not affect the ionization history at z > 1600, 
while at z ~ 1000 one has ANJN^ ~ -\AP/P. Also one can see 
that at low redshifts, changes of the escape probability are not 
propagating very much to the ionization history, resulting only 

in ANe/Ne 0.054 Af/P at z ~ 600. Note that the correction 

of Ne; in both cases is much smaller than the one of the escape 
probability. This is because the 2s- Is two-photon decay channel 
already contributes slightly more to the effective recombination 
rate. 

In Fig. [12] we give the coiTection of the number density 
of free electrons as a function of redshift. The time-dependent 
correction of the escape probability leads to a ~ 1.6 - 1.8% 
change of ANe/Ne in the redshift range 800 < z < 1200. 
This chang e is practically twice as lar ge as the effect due to 
line recoil (iGrachev & Dubrovichll2008l) . so that the sum of the 
time-dependent correction and the recoil coiTection is still dom- 
inated by the former contribution, leaving ANe/Ne ~ 0.94% at 
z ~ 1 100. This will be important for the computation of the CMB 
power spectra, where at large / in the case of TT one expects a 
ACi/Ci ~ 1%, and about 2 times more for EE. 
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5. Discussion and conclusions 

5.1. Main results related to the cosmological ionization 
history and the CMB power spectra 

In this paper we investigated the validity of the Sobolev approx- 
imation for the Lyman a escape probability during hydrogen re- 
combination. We separate absorption and emission of Lyman a 
photons from resonant scattering events, including the fact that 
processes leading to full redistribution of photons over the Voigt 
profile occur with much lower (~ lO"-' - 10"'* times) probability 
than resonant scatterings. We have shown that within the stan- 
dard formulation the rapid changes in the ionization degree dur- 
ing recombination lead to significant departures of the photon 
distribution from the quasi-stationary solution. We took these 
coiTections into account analytically, assuming that the photon 
redistribution process over frequency during a scattering is co- 
herent in the lab frame. 

Although one does expect some additional modifications 
when accounting for partial frequency red istribution, our com- 
putations show (Chlu ba & SunvaevI 12009 ) that the additional 
correction will be do minated by the influence of line recoil, that 
has been addressed in iGrachev & DubrovichI (|2008|) . However, 
the time-dependent correction that is considered here turns out 
to be significantly larger, so that we focused on this only. A 
more complete consideration of this problem is in preparation 
dChluba & Sunyaev 2009) . 

Here we found that the time-dependent corrections to the ef- 
fective Lyman a escape probability result in a ~ 1.6 - 1.8% 
change of ANe/Ne in the redshift range 800 < z < 1200 (See 
Sect. 13.41 andl4lfor more detail). These corrections are impor- 
tant for the Thomson visibility function and in computations of 
CMB power spectra, where at large (I ~ 1000 - 3000) multi- 
poles / in the case of TT one expects modifications of the order 
of ACi/Ci ~ 1%, and about 2 times more for EE. However, note 
that we also expect additional changes when formulating the 
problem more rigorously in the two- or multi-photon approach 
(see discussion in Sect. 15.21 ). 

The main reason for the corrections discussed here are (i) 
time-dependent changes in the absorption optical depth, and (ii) 
changes in the net emission rate due to the time-dependence of 
cosmological recombination. The coiTection due to case (i) is 
especially important for contributions coming from the distant 
wings of the Lyman a line, where emitted photons can travel, 
scatter, and redshift over a very long time before getting reab- 
sorbed. For the 3 shell hydrogen atom the associated correction 

to the escape probability is AP/P 1 .2% at z ~ 1 100, while it 

results in an increase of AP/P — h2.8% at z ~ 1440 (for more 
details see Fig.|5]l. 

The correction related to the time dependence of the net 
emission rate is slightly larger, leading to a total decrease of the 
effective escape probability at z < 1400, that for the 3 shell atom 
reaches AP/P ~ -5.8% at z ~ 1140, while it results in an in- 
crease of AP/P ~ +3.7% at z ~ 1510 (for more details see 
Fig. nil. Here a significant contribution is due to departures of 
the very distant (~ 10^ - lO'* Doppler widths) wing spectrum 
from the quasi-stationary solution (e.g. see Sect. l3.4.5] l. 

We also showed that in particular at low redshifts (z < 
800 - 900) this correction, owing to a self-feedback process, 
strongly depends on the number of shells that have been taken 
into account for the computation (see Fig. This is because 
the effective emission rate depends on the solution for the popu- 
lations of the excited levels, so that it is important to include at 
least 4-5 shells into the computations. However, this aspect of 
the solution appears to be due to the incompleteness in the for- 



mulation of the problem, so that at these redshift the conclusions 
should change when using a two- or multi-photon description 
(see discussion in Sect. 15.2b . Also it is important to mention that 
for the corrections in the CMB power spectra, this should not 
affect the results very much. 

5.2. Apparent problems with the standard formulation 

Our analysis shows that under the extreme physical conditions 
valid in the hot Universe (extremely low plasma density in the 
presence of the intense CMB radiation field), the standard for- 
mulation of the Lyman a transfer problem leads to several ap- 
parently unphysical results. First, we would like to point out that 
all our computations and estimates clearly show how important 
(at the percent level accuracy) the distant wings of the lines are 
for the value the escape probability or mean intensity supporting 
the 2p level (e.g. see Sect l3.3.4l and l3.4.5] l. However, in the stan- 
dard approach variations of the blackbody and also any power- 
law variations in v are usually neglected in the formulation of 
the transfer problem and analytic computation, an approxima- 
tion that is certainly questionable when going to |Av/v| > 1%, 
or ~ 10^ Doppler width. For example, as we mentioned in 
Sect. 12.2.41 this approximation leads to a small non-conservation 
of a blackbody spectrum at large distances from the line center, 
an aspect that simply follows from exact application of the de- 
tailed balance principle, leading to a thermodynamic coiTection 

factor / - ^ in a two-photon formulation of the 

proble m this factor automatically appears dChluba & SunvaevI 
'2008a'). 

Also, the emission of photons according to the standard 
Voigt profile in principle allows the production of photons un- 
til arbitrarily large distances on the blue side of the Lyman a 
resonance. Without introducing some high frequency cut-off, 
in the cosmological recombination problem these photons will 
lead to some unphysical self-feedback at low redshift (e.g. see 
Sect. 13.4.1b . which is also present in our current solution, but at 
times that are not so important for the CMB power spectra. We 
expect that both problems can be resolved when using a two- or 
multi-photon formulation, in which detailed balance is applied 
self-con sistently, and where the line profiles are naturally bound 
(e.g. see lChluba & Sunvaevll2008bl) due to energy conservation. 

Focusing on the Sobolev approximation (quasi-stationarity 
of the spectrum and complete redistribution), several unphysical 
aspects also appear. These are again due to the unique properties 
of our Universe, where there are hardly any collision and the ex- 
pansion rate is so low that the Sobolev optical depth ts reached 
values of ~ 10^-10** during recombination. As explained in 
Sect. |3.2| this leads to the case that the variations of the photon 
distribution that are important for the mean intensity support- 
ing the 2p level during cosmological recombination occur at dis- 
tances of ~ 10^ - 10'^ Doppler widths from the resonance. This 
is far beyond the Lyman p line or even the ionization energy of 
the hydrogen atom. 

It is also possible to compute the present day Lyman a spec- 
tral distortion in the time-dependent approach, using the solu- 
tion ( |26] |. for which it was assumed that every line scattering 
leads to a complete redistribution of photons over frequency. We 
checked that in this case one would obtain a Lyman a line pro- 
file that is very different from the one computed in the usua l 
6 function approximation (e.g. see Rubino-Ma rtm et al.ll2006l) . 
One reason for this is that the effective frequency beyond which 
the photon distribution is no longer affected by the Lyman a res- 
onance is very far on the red side of the line center (at xd 10^ 
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to -10"^, or |Av/v| ~ 1% - 10%). This aspect of our computa- 
tions also suggests that in principle it should be possible to con- 
strain the type of redistribution that is at work during hydrogen 
recombination by looking at the exact position and shape of the 
residual Lyman a distortion in the CMB. In all cases, no line 
scattering, complete redistribution, and partial redistribution the 
Lyman a distortion will look different. 

We conclude that for the conditions during cosmological re- 
combination, complete redistribution for a line scattering event 
in the standard formulation is not an appropriate redistribution 
process, and leads to rather unphysical results. With the ap- 
proximation of coherent scattering in the rest-frame, some of 
the unphysical aspects of the solution disappear However, as 
mentioned above, at low redshifts we obtain a large feedback 
of Lyman a photons initially released at high redshifts (e.g. see 
Sect. 13.4.1b . These problems can be resolved using a two- or 
multi-photon formulation. 



apply the Taylor expansion 



^wings 9 



15-20a2 105(1 -2a2) 



2x^ 



(A.l) 



We checked that the Voigt function is represented with relative 
accuracy better than 10"^ in the whole range of frequencies and 
redshifts. Using Eq. lA.ll on the red side of the resonance one 
can approximate the integral - J^^ tp{x') Ax' by: 



X 



wmgs 



1 + 



3 - 2fl2 3 - 4fl2 



15(1 -2fl2) 



6x1 



(A.2) 



as long as x < -30. Since a ~ 10"^ - 10"^, this shows that in the 
distant wings only a very small fraction of photons is emitted. 
Using the symmetry of the Voigt profile one finds x(x) = 1 - 
Xi-x), such that Eq. ( IA.2b is also applicable for x > 30. 



5.3. Future prospects 

In spite of all the complications we expect that to lowest order 
one can take the time-dependent correction during cosmolog- 
ical recombination into account using the solution (|49] |. or in 
the full two-photon description using a time-dependent solution 
in the no-scattering approximation. Since all additional correc- 
tions will (also) be small, one can then compute each other pro- 
cess more or less separately. This should also be possible for the 
case of helium recombination, but here the reabsorption of pho- 
tons by the small fraction of neutral hydrogen atom s present at 
that time will be much more important ( S witzer & H iratal2008at 
iKholupenko et al.ll2007l:lRubino-Martm et al.,.2008.) . 

In order to include the final correction into the computa- 
tions of the CMB power spectra it will be necessary to develop a 
fast scheme for the evaluat ion of the ionization history. For this 
purpose Fen dt et al.l (l2008h recently proposed a new approach 
called RiccP^ which uses multi-dimensional polynomial regres- 
sion to accurately represent the dependence of the free electron 
fraction on redshift and the cosmological parameters. Here one 
first has to produce a grid of models using a given full recom- 
bination code, for which each run may take several hours or up 
to days. However, the time-consuming part of the computation 
is restricted to the training of Rico, while afterwards each call 
only takes a small fraction of a second. This approach should 
allow one to propagate all the corrections in the ionization his- 
tory that are included in the full recombination code to the CMB 
power spectra, without using any fudge-factors, like in Recfast 
dSea^er et al] 1 19991: IWong et al.li2008) . In the future, we plan 
to provide an updated training set for Rico, including the time- 
dependent corrections discussed here. This should also make it 
easier for other groups to cross-validate our results. 



Appendix A: Computational details 

A. 1 . Computations of 0(v) 

The evaluation of the Voigt profile, Eq.(fT3]l, is usually rather 
time-consuming. However, convenient approximations for 4>{xd) 
can be given in the very distant wings and also close to the center 
of the line. For \xq\ < 30 we use the approximati on based on the 
Dawson integral up to sixth order as described in lMihalai (Il978l 
Sect. 9.2, p. 279). In the distant wings of the line (|xdI > 30) we 
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